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We develop a technique based on density functional theory and the projector augmented wave 
method in order to obtain the x-ray absorption cross section at a general edge, both in the electric 
dipole and quadrupole approximations. The method is a generalization of Taillefumier et aZ..— We 
apply the method to the calculation of the Cu L,2,3 edges in fee copper and cuprite (CU2O), and to 
the S L,2,3 edges in molybdenite (2H-M0S2). The role of core-hole effects, modeled in a supercell 
approach, as well as the decomposition of the spectrum into different angular momentum channels 
are studied in detail. In copper we find that the best agreement with experimental data is obtained 
when core-hole effects are neglected. On the contrary, core-hole effects need to be included both in 
CU2O and 2H-M0S2. Finally we show that a non-negligible component of S L,2,3 edges in 2H-M0S2 
involves transition to states with s character at all energy scales. The inclusion of this angular 
momentum channel is mandatory to correctly describe the angular dependence of the measured 
spectra. We believe that transitions to s character states are quantitatively significant at the L2,3 
edges of third row elements from Al to Ar. 
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I. INTRODUCTION 

X-ray absorption near edges structure (XANES) spec- 
troscopy is a very powerful tool for condensed matter 
studies. Its chemical and spatial selectivity allow solving 
crystallographic structures and investigating the proper- 
ties of materials, such as magnetism or orbital hybridiza- 
tion. Unambiguous interpretation of x-ray spectra re- 
quires theoretical insight, in order to determine the ori- 
gin of spectral features and to understand the effect of 
the core-hole on the measured spectrum. 

Several theoretical approaches, cither in real spaced 
or based on periodic boundary conditions^— were proven 
successful to describe XANES spectra in molecules and 
solids at the density functional theory (DFT) level. 
All-electron based implementations^ the most accurate 
within DFT, suffer from significantly larger computa- 
tional costs than those required by pseudopotential- 
based methods.— ~— In this respect, the advantage of 
pscudopotcntial-bascd methods is that they allow to 
tackle larger systems, e.g. composed of several hundreds 
of atoms. In pscudopotcntial-based methods the atomic 
core-potential is replaced by a fictitious one, smoother 
and less singular, hence requiring a smaller representa- 
tion in terms of plane waves and consequently less com- 
putation time. As a consequence, in order to describe 
properties that depend on core-states, like XANES, one 
needs to reconstruct the all-electron wavefunction of the 
empty valence states. This is achieved in the framework 
of the projector augmented wave (PAW) method^ 

The PAW approach for the calculation of x-ray absorp- 
tion spectra has been implemented in some freely avail- 
able (for K-edge only)£ first-principles codes, commer- 
cial code distributions 7 -^ and in other non-distributed 
software. 8 ' 11 ' 12 The PAW method has been proven suc- 



cessful in interpreting experimental dat a 8 ' 11 ' 13 " —, thanks 
to its ability to perform the spectral analysis of the cal- 
culated XANES. Most of the applications of the exist- 
ing methods involve K-edge calculations, where usually 
core-hole is reliably modeled in a supercell approach and 
electronic excitations are well described by DFT. 

PAW studies at L*2 : 3 edges are not numerous . 7 ' 22 ' 23 
One of the main reasons is that, contrary to K edges, 
many-electron effects are common at L.2,3 edges. For in- 
stance, in the case of L2.3 edges of transition elements and 
M4 : 5 edges of rare earths, excitations cannot no longer 
be treated at the DFT level, due to the breakdown of 
the single-particle pictured Consequently, more sophis- 
ticated methods based on time dependent DFT^i 2 ^ mul- 
tichannel multiple scattering theor y 27 ' 28 or the Bcthe- 
Salpctcr equation 2 ^ need to be employed. However, none 
of the aforementioned first principles methods can rigor- 
ously deal with open shells. In this particular case, the 
ligand field multiplet theory 3 ^ (LFM) can successfully in- 
terpret the spectra, albeit not from first principles. 

When the states probed by the spectroscopy are delo- 
calizcd - K and Li edges in general (with the exception 
of very light materials), L2.3 and M^ edges of elements 
belonging to the Ad and 5d series - x-ray absorption spec- 
tra can usually be interpreted in the DFT framework. In 
these cases, the final states, mainly of p (K and Li) or d 
character (1^3 and M2,3 edges) show limited correlation 
effects and, as a consequence, the structures in a XANES 
spectrum can be associated to the Kohn-Sham states. 

In this work we focus on the description of 1^3 edges 
within DFT and PAW frameworks. The previous work 
in Ref. [0 is generalized to the case of a general edge, for 
a spin-orbit split core-level. In a first step, we validate 
our method by calculating Cu L 2 ,3 edges of metallic Cu. 
We study the convergence of the spectra with respect to 



the number of projectors considered in the reconstruc- 
tion step, i.e. the completeness of the projector basis. It 
is generally agreed^ that at K edges two linearly inde- 
pendent projectors per angular momentum channel are 
sufficient as to insure an accurate reconstruction of the 
Is state. Albeit a pre-requisitc for PAW applications, 
other edges, and L2.3 in particular, have not benefit from 
similar studies. 

Next, we switch to the more challenging cases of the Cu 
L2.3 edges in cuprite (CU2O) and S L 2j 3 edges in molyb- 
denite (2H-M0S2). For both compounds, we perform a 
detailed analysis of the electronic origin of the spectral 
features and we assess the effect of the 2p core. Concern- 
ing CU2O, we investigate how the empty states probed 
by the spectroscopy are influenced by the inclusion of 
the Hubbard U term. In the case of the hexagonal 2H- 
M0S2 we interpret the effect of the polarization of the 
beam, i.e. the x-ray natural linear dichroism (XNLD), 
by decomposing the spectra into various angular moment 
channels. 

The article is organized as follows. In section HI1 we re- 
call formal aspects of the x-ray cross section calculation in 
the PAW formalism. In section ITO1 we give the technical 
details of the calculations. Next, in section lTVl we present 
results for the three aforementioned compounds. Finally, 
in the appendix we include the full analytical expression 
of the electric dipole (El) and electric quadrupole (E2) 
matrix elements for a general edge. 



II. X RAY ABSORPTION CROSS SECTION 

In the single particle approximation, the x-ray absorp- 
tion cross section reads ^ 

a(u) = An 2 ahwY^ I (ip f \d\ipi) | 2 S (Hu - E f + E % ) (1) 

where ipi is the one-particle wavefunction of the core level 
and ipf is the all-electron one-particle (empty) final state. 
The corresponding electronic energies are Ei and Ef, hui 
being the photon energy and a is the fine structure con- 
stant. The transition operator is taken as: 

6 = V+Q 

where T> = i ■ r and Q — |(e • r)(k ■ r) are the electric 
dipole and quadrupole operators, respectively. The quan- 
tities, e and k are the polarization vector and the wave- 
vector of the incident x-ray beam, respectively, and r is 
the position coordinate of the electron. In this work we 
only deal with the pure dipolar absorption term (E1E1). 
The quadrupolar contribution (E2E2) is generally negli- 
gible at the 1^2,3 edges of transition elements, but may be- 
come important in rare-earths compounds. E1E2 cross- 
terms appear in optically active materials. 

In the next paragraphs we briefly introduce the PAW 
method applied to the calculation of the x-ray absorption 
cross section. We closely follow Refs. [l| andP 



In a pseudopotcntial calculation the pscudo- 
wavefunction of the crystal \ipf) is obtained at the 
end of the self-consistent field run. In order to get 
the all-electron wavefunctions \tpf) needed in Eq. Q] 
all electron reconstruction^ needs to be performed. As 
shown in Refs. [0, and |32| the cross section may be 
written as: 



c(w) = 4-7T atnj y 
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In Eq. [3l |</>r , p ) {\(f>R , P }) are the all-electron (pseudo) 
partial waves centered on the absorbing atom. The quan- 
tities (pr 0)P | are a complete set of projector functions, 
labeled by the index p. The following conditions are sat- 
isfied: 



^r,p( r ) = fe,pM outside Or 






(Pr,p <hu,p l ) = <*RR' <W inside Or 



(4) 
(5) 



The quantity Or, is the so-called augmentation (or core) 
region centered on atomic sites R. In our case, the all- 
electron partial waves are chosen to be the solutions of 
the Schrodingcr equation for the isolated atom. 

It is worthwhile recalling that Eq. [2] holds under the 
assumption that the initial state is localized on the ab- 
sorbing atom, so that the overlap between the state 0\ipi) 
and \<Pr) can be neglected if R 7^ Ro- Furthermore, while 
Eq. [3]is in principle valid for a complete set of projectors, 
it converges after a few terms. 

In the appendix section [A] we derive the expression of 

5 R ,p ^ ipi) f° r a general 



the transition matrix element < </>r 0iP O 

edge. 

The single particle interpretation of x-ray absorption 
is no longer pertinent in the presence of a strong in- 
teraction between the photoelectron and the core hole. 
For the E1E1 (electric dipole) channel, this is the case 
when the spin-orbit splitting of the initial state is small 
and the adjacent edges overlap. A notable example are 
the L-2,3 edges of transition elements 2 ^ or M^ edges of 
rare earths, no longer obeying the statistical branching 
ratioi 2 ^ At K edges, neglecting the interaction between 
the photoelectron and the core hole leads to the over- 
estimation of the energy position of the E2E2 (electric 
quadrupole) spectral structured Often, the angular de- 
pendence of the E2E2 structure is inappropriately de- 
scribed in the single particle pictured 



III. TECHNICAL DETAILS 

The calculation of a general edge was implemented in 
the XSpectra codo£ belonging to the Quantum Espresso 



distribution^ The method was applied to the calculation 
of x-ray absorption spectra of fee Cu (L 2j3 edges of Cu), 
Cu 2 (L 2 ,3 edges of Cu) and 2H-MoS 2 (L 2 ,3 edges of S). 

The experimental lattice parameters arc used for all 
three systems: a = 3.600 A for fee Cu^ 4 - (spaccgroup 
Fm3m), a = 4.268 A for Cu 2 Q^ (spacegroup Pn3) and, 
finally, a = 3.161 A and c = 12.295 A for 2H-MoS2^ 
(spacegroup P63/mmc). 

For Cu, Mo and O, we used ultrasoft 
pscudopotcntiala^ including scmicore states for Cu 
and Mo. We use a Troullier-Martin^ norm-conserving 
pseudopotential for S. Three, respectively two projectors 
for the I = 2 channel, as well as one projector for the 
/ = channel were needed for the proper reconstruction 
of core states for the absorbing Cu and S atoms. 

Cu and Cu 2 are described within the generalized 
gradient approximation (GGA) in the Perdew-Burke- 
Ernzerhof (PBE) paramctrization^ In the case of 2H- 
MoS 2 we used both the local density approximation 
(LDA) via the Pcrdcw-Zunger (PZ) functional 4 ^ and 
PBE, without any detectable difference. 

The kinetic energy cutoff of the wavefunction expan- 
sion is chosen as 30 Ry in Cu, 60 Ry in Cu 2 0, and 50 Ry 
in 2H-MoS 2 . The charge density cutoff is set ten times 
larger then the kinetic energy one. 

In the case of metallic Cu, integration over electron- 
momentum is performed by using a uniform 20 x 20 x 
20 grid. In the case of insulating Cu 2 and 2H-MoS 2 we 
use 20 x 20 x 20 and 10 x 10 x 5 electron-momentum 
grids, respectively. 

The core-hole is described within a static approxima- 
tion. A pseudopotential with a core-hole is generated 
for the absorbing atom. Furthermore, a supercell is used 
to minimize the interaction between the absorbing atom 
and its periodic images. We use a 3 x 3 x 3 supercell 
for Cu, a 2 x 2 x 2 one for Cu 2 and finally a 2 x 2 x 1 
supercell for 2H-MoS 2 . 

In the XANES calculations, electronic integration was 
performed on a 30 x 30 x 30 grid for Cu and on a 20 x 
20 x 10 grid for 2H-MoS 2 . In the case of Cu 2 we use 
the same grid as in the self-consistent calculation. These 
grids refer to the unit cell when no core-hole is involved. 
For supcrcells the grids are rescaled, due to the smaller 
Brillouin zone. 

The spin-orbit coupling on the final (valence) states 
probed by the spectroscopy is neglected. Although full 
relativistic calculations arc needed to describe heavy 
transition elements of the late Ad and the bd series, 
our approximation is reasonable for the applications pre- 
sented in this paper. In particular, the spin-orbit cou- 
pling acting on the 3d states of transition elements sub- 
jected to crystal field is negligible! 41 ' 42 However, the de- 
scription of the spin-orbit coupling for the core states is 
mandatory in order to distinguish between related edges 
as the L 2 (transitions from 2pi/ 2 ) and L3 (2p 3 / 2 ), M 2 
and M3 or M4 and M5. In the appendix we show the 
transition matrix elements calculated accordingly. 

The radial parts of the 2pi/ 2 and 2p 3 / 2 core wavefunc- 
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FIG. 1. The raw spectrum is convoluted with a lorentzian 
of acrtangent-like, energy-dependent width. Ef is the cutoff 
energy, Thole is the spectral width at low energy, above Ef 
and Y max is the value at high (infinite) energy. The inflection 
point is situated at E c t r . 



tions are taken as identical, a very good approximation 
in the applications treated in the present paper, as it 
can be explicitly verified from a DFT relativistic all elec- 
tron calculation on the isolated atom. We account for 
the spin-orbit splitting of the 2p states by shifting the 
calculated spectra by the value of the energy separation 
between the L3 and the L 2 edges. This value can be 
either taken from a DFT relativistic all-electron calcula- 
tion on the isolated atom or from experiments. The two 
choices are equivalent when the core-hole can be treated 
at a DFT level. On the contrary, strong Coulomb inter- 
action between the 2p hole and the valence d electrons 
may lead to an energy separation between the L3 and the 
L 2 which is different from the value of the core level spin- 
orbit splitting^ In the absence of spin-orbit coupling on 
the d states, the L 2 and L3 edges have similar structures, 
while the intensities are related by a factor of 2 (the sta- 
tistical branching ratio) , according to the manifold of the 
initial states. 

The occupied states are eliminated from the XANES 
spectra by introducing a cutoff, which appears implicitly 
in Eq. [1] in the summation on the empty states /. The 
cutoff is chosen at the calculated Fermi level or at the 
lowest unoccupied band, depending on the metallic char- 
acter of the compound. In the presence of a core-hole, 
the self-consistent calculation is performed on a charged 
cluster. This means that the energy cutoff is calculated 
by assuming that the photoelectron does not participate 
to the screening of the core-hole. 

The XANES spectrum being calculated from station- 
ary Kohn-Sham states, a broadening must be applied in 
order to account for the finite life-time of the electronic 
levels involved in the transition: 
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The Fermi level or the lowest unoccupied band energy 
is labeled Ep. Subsequently, the raw spectrum a(uj) is 
convoluted with a lorentzian having an energy dependent 
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FIG. 2. (color online) L2,3 absorption edges in metallic Cu. 
Experimental data— (dots) versus calculations: without core- 
hole (thick solid) and with half core-hole (thin solid). The 
spectra are shifted vertically for convenience. 



FIG. 3. Two distinct models for the convolution of spectra. 
Lower panel: energy independent (dashes) and arctangent- 
like Lorentzian- width (thin solid). Upper panel: the corre- 
sponding convoluted spectra for Cu L2.3 absorption edges: at 
constant width (dashes) and arctangent-like (thick solid). 
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The first term Thole is the core-level width, taken as 
energy-independent since screening is treated statically. 
In practice we use the standard, edge-dependent, tab- 
ulated values4£ The second term 7(f2) is the spectral 
width due to the final state. In general, the spectral 
structures at the L 2 edge are broader than the ones at 
L3 due to the presence of an extra desexcitation channel, 
the Coster-Kronig effect^ Hence, distinct convolution 
widths are needed at the two edges in order to repro- 
duce the spectral shape. Furthermore, to account for 
inelastic scattering events (e.g. plasmons), a Lorentzian 
with a more sophisticated energy dependence is some- 
times preferred to a step-like function. We choose the 
convolution with an arctangent-like function, an empiri- 
cal model close to the Scah-Dcnch formalism^ 
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where e = (Q — Ep)/(E ctr — Ef). The energy dependent 
broadening is bound between Thole and T max , the values 
at low and respectively, high (infinite positive) energies. 
The inflection point of the arctangent is situated at E ctr . 
Figure Q] illustrates the generic form of r/(0) and the 
significance of the above parameters. In this work, the 
following values of Thole, T max and E ctr ~ Ep have been 
used: 0.35, 4, 6 eV for Cu; 0.22, 4, 9 eV for Cu 2 and 
0.3, 7, 6.5 eV for 2H-MoS 2 . 



IV. RESULTS AND DISCUSSION 
A. Cu L2,3 edges in copper 

To validate our implementation we first calculate L2.3 
edges for fee Cu. In fee Cu, the 3d states of the absorb- 
ing atom are almost entirely occupied (atomic Cu is in 
nominal 3d 10 ). XANES at L 2 3 edges probes the first 
empty d states just above the main 3d occupied band. 
These states contribute very weakly to the overall DOS 
and in a narrow energy region. For this reason, and in 
agreement with the result in Ref. 1471 . we find that an 
extremely dense k-point sampling (electron-momentum 
grid) is needed to converge the x-ray spectra (see the 
technical details section). 

In Fig. [2] the calculated x-ray absorption (XAS) spec- 
tra are compared to experimental data4£ The effect of 
the core-hole at the Cu L 2]3 edges in fee Cu was thor- 
oughly tested. We find that the best agreement between 
theory and experiment is achieved by neglecting core-hole 
effects. The global intensity mismatch at the L 2 edge (the 
calculated intensities are less than the measured one) is 
most probably due to the background subtraction in the 
experimental spectrum. The inclusion of the full core- 
hole (not shown) shifts the weight of the d states at sev- 
eral eV below the Fermi level and significantly decreases 
their contribution above the Fermi level. Even with half 
core-hole, we find that the agreement with the exper- 
iment is worse than in the case without core-hole, with 
the first peak at the L3 being less intense than the second 
one (see Fig. [2]). 

Our analysis is in partial agreement with the previ- 
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FIG. 4. Dependence of the calculated Cu L2 edges on the 
core- level reconstruction: one (dashes), two (thin solid) and 
three (thick solid) projectors for the I — 2 channel. 



ous interpretation of XANES data on fee Cu in Ref. |48|, 
based on the Korringa-Kohn-Rostoker method of band 
structure calculation. The only notable difference be- 
tween the two calculations is that Ref. |48| produces full 
relativistic results. In both calculations, the core-hole 
was not included explicitly. While the positions of the 
spectral structures agree in the two calculations, it is no 
longer the case for the intensity of the leading L3 peak. 
More precisely, the ratio between the leading peak and 
the second L3 structure in the spectrum of Ref. |48| is 1.5 
larger than in our calculations. 

The reason for this disagreement is very likely related 
to the procedure adopted to eliminate the occupied states 
from the calculation. We find that the intensity of the 
first peak in the spectrum strongly depends on this fac- 
tor, which is often a problem in metallic systems, due to 
the absence of the gap in the electronic spectrum. As 
such, the disagreement between the two calculations is 
not due to fundamental methodological differences in the 
calculation. 

Another difference with Ref. |48| concerns the estimate 
of the magnitude of transitions from 2p to s states. In 
Ref. |48| it was estimated to be 5% of the total intensity, 
whereas they are completely negligible according to our 
results. The discrepancy may come from treatment of 
relativistic effects. 

Note that the analysis of electron energy loss near edge 
structure spectroscopy (ELNES) data ; 47 i 49 carried out 
with a similar formalism to the one used in XANES, leads 
to the conclusion that a partial half core-hole screening 
is necessary to explain the experimental data. This ap- 
parent disagreement is mostly due to the differences in 
the experimental XANES and ELNES spectra, and not 
to our calculation, which actually agrees fairly well with 
the ones featured in Refs. |47| and |49|. Indeed, the first 
peak L3 peak in ELNES spectr a 47 i 49 is smaller than the 



second one, opposite to XANES experimental dat a 46 i 48 
(see, for instance, Fig. [5]). 

The fact that we achieve a very good description of 
the experimental data by completely neglecting the core- 
hole should be interpreted carefully. While one might be 
tempted to say that in metallic systems the core-hole is 
almost entirely screened, we recall that at K edges this is 
not the casei&££ It has already been pointed oul^i that it 
is unlikely to predict whether a XANES calculation with 
core-hole will agree better to the experiment than the one 
without hole, due to the inadequacy of the description of 
the core-hole in the single particle approximation. Nev- 
ertheless, one can say that the 2p hole is less profound 
than the Is one and consequently less screened. There- 
fore partial screening is more likely to be needed for the 
description of L 2 ,3 edges than at K edges. 

Fig. [3] features the effect of the convolution procedure 
on the shape of L2,3 edge of Cu (see Eq. [7|). The in- 
clusion of an energy dependent convolution improves the 
agreement with the experiment at the L2 edge. Indeed, 
a Lorentzian of constant width {Thole) overestimates the 
resolution at high energies, whereas the arctangent-like 
energy-dependent Lorentzian can reproduce the correct 
broadening. 

We recall that the sum in Eq. [3] formally involves an 
infinite number of terms. However, experience demon- 
strated that at K edges, two linearly independent pro- 
jectors for the I = 1 channel are enough in order to con- 
verge the E1E1 calculated XANES within the first « 50 
eV above the edge£ The description of the E2E2 contri- 
bution sometimes requires three projectors for the I = 2 
channel, particularly when the energy of the first projec- 
tor is set at the one of the semicore states. It is then 
not obvious that two I = 2 projectors are sufficient for 
calculations of L2,3 edges, where the sampled states have 
mostly d character. 

Fig. 2] shows the convergence of the L2 spectrum with 
respect to the number of projectors with d character con- 
sidered. In particular, we set the first projector at the 
energy of the DFT all-electron 3d states of atomic Cu in 
neutral configuration. The second and third projector are 
at 2.5 Ry and 10.0 Ry higher energies. This choice leads 
to linearly independent projectors. As shown in Fig. |4l 
the use of a single projector underestimate the XANES 
cross section already at energies just above the main L2 
edge. The use of additional projectors increases the cal- 
culated L2 XANES at higher energies. Three projectors 
are by far sufficient to describe the full L2 spectrum 50 
eV above the edge. 

In Ref. [D the authors recall that Eq. [5] is not nec- 
essarily exact for the 2p core level. Indeed, should the 
2p wavefunction no longer be confined within the aug- 
mentation region, Eq. [2] is to be replaced with its more 
general form (see, for instance, Eq. 5 in Ref. |l|). How- 
ever, we find that in the case of Cu both Is and the 2p 
core wavefunctions are well localized within the augmen- 
tation region, hence the reconstruction according to Eq. 
[5] is undoubtedly correct, provided that (i) a sufficient 
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FIG. 5. (color online): electronic structure of CU2O: PBE 
(crosses) and PBE+U (circles) with U — 8 eV. In the right 
panel, the electronic DOS is plotted: PBE (thin solid) and 
PBE+U (thick solid). Positive energies are associated to 
empty states. 



number of projectors is included and (ii) the radius of 
the augmentation region Qr is large enough to confine 
the 2p radial function. 



B. Cu L,2,3 edges in cuprite (CU2O) 

Having validated our implementation on fee Cu, we 
switch to the more complex case of CU2O. CU2O is a non 
magnetic insulator where copper is nominally in a Cu +1 , 
3d 10 valence state. Density functional theory (DFT) suc- 
ceeds in capturing the insulating ground state of CU2O, 
albeit with a smaller gap than the one detected in exper- 
iments. Our PBE calculations predict a 0.36 cV direct 
gap, in agreement with previous calculations ^2r— to be 
compared to the much larger experimental value, rang- 
ing from 2.02 to 2.4 ± 0.3 eV (see Ref. [55| and references 
therein) . 

A way to improve the description of the gap is to in- 
clude the Hubbard U in the calculations, via the DFT+U 
scheme. While this has been achieved in Ref. [5(| the 
question whether it improves the description of the empty 
states seen in the x-ray absorption spectroscopy is still 
open. 

Fig. [S] features the band structure over selected high- 
symmetry directions in the Brillouin zone, for U = 8 
cV. One can see a small opening of the gap, from 0.36 in 
PBE to 0.43 eV in PBE+U. The reason for the small gap 
opening has been explained in Ref. l56l where the authors 
acknowledge that both DFT and DFT+U overestimate 
the amount of s-d hybridization. 

While the gap increase is very small, the empty states 
are substantially shifted to higher energy after inclusion 



of the U term (Fig. [5]). This will significantly affect the 
shape of the calculated XANES spectra, which can be 
safely interpreted within the single particle picture, given 
the (nominally) closed shell configuration. Furthermore, 
by direct comparison with the experiment one can assess 
whether the DFT+U approximation is appropriate to the 
description of CU2O. 

In order to see the qualitative effect of the Hubbard 
term, we choose the value of U = 8 eV that is (i) a 
realistic value for Cu based oxides&iS and (ii) allows a 
substantial shift of spectral features to judge the effect 
of U in this system. 

CU2O crystallizes in a cubic structure with four equiv- 
alent Cu atoms per unit cell. To calculate the absorp- 
tion signal, one should consider the contribution of all 
Cu atoms (one calculation for each absorber). However, 
given that the four Cu atoms in the unit cell are equiv- 
alent, the problem can be simplified by making use of 
symmetry considerations. Let the atomic E1E1 absorp- 
tion tensor d defined by 
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where a, j3 = x,y,z and i = 1,2,3,4 runs over the 4 
equivalent atoms. The matrix d^> is hermitian and obeys 
the point group symmetry of the specific atom site. In the 
case of the Cu atom located at (j \ j) in reduced crystal 
coordinates (first conventional origin choice for Pn3), the 
corresponding tensor is (Se point group symmetry): 



d« 



with a 7^ b real. The remainder d«) with j = 2,3,4 can 
be obtained by: 




d« 



74 d« n l3 



where IZij is the 3x3 cartesian rotation matrix associ- 
ated to the symmetry operation relating the prototypical 
atom i to its equivalent j. In our case the (j \ j) Cu 
atom is related by C3 rotations to the other Cu atoms. 
The crystal absorption tensor V can be obtained as the 
average over all the atomic tensors: 



'Dap = 



\t 



+) 

'ail 



In matrix form: 



a 

I a 

0a 



(8) 



(9) 



Note that V is isotropic, which was to be expected for the 
absorption tensor of a crystal described by a cubic space- 
group. On the contrary, the atomic absorption tensors 
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FIG. 6. (color online) Calculated XANES spectra (solid) at 
the Cu L<2,3 edges in CU2O and comparison with experimental 
data from Ref. |57| (dots): PBE with no core- hole (top) and 
with half core-hole (bottom). 
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FIG. 7. (color online) Calculated XANES spectra (solid) at 
the Cu Ij2,3 ed ges in CU2O and comparison with experimental 
data from Ref. [57] (dots): PBE+U with no core-hole (top) and 
with half core-hole (bottom) 



S l > are not isotropic, a consequence of the Sq point group 
symmetry. Hence we calculate the total E1E1 absorption 
signal per Cu atom as the absorption of the prototypical 
(2 4 4) *-' u a t° m m the (100) direction. 

Similar to the case of Cu, we performed XANES calcu- 
lations with three distinct screening schemes: with core- 
hole, with half core-hole and without core-hole. The re- 
sults are shown, with or without the inclusion of the Hub- 
bard U, in Figs. (U = 8 eV) and© (U = eV) respec- 
tively. Calculations were normalized to the edge-jump 
at high energies and aligned to match the experimental 
peak at 934 eV. Similarly to the case of Cu, we find that 
the inclusion of a full core-hole (not shown) shifts the d 
states completely below the Fermi level and hence yields 
unphysical results. In Fig. [7] wc show that the inclu- 
sion of the Hubbard U term in the description of the d 
states of Cu in CU2O gives a very poor agreement with 



the experimental data, independently of the screening, 
contrary to what has been observed in other oxides^ 
We thus bring supplemental evidence to the finding in 
Ref. [56| according to which GGA+U does not contain 
the proper physics to describe Cu 2 0. 

The inclusion of a half core-hole is enough in order to 
explain the measured XANES at the L 2i 3 edges of Cu. 
In Fig. |51 we show that calculations without core-hole 
at U = eV overestimate the L3 structure at 931 eV. 
On the contrary, a good agreement with the experiments 
is achieved with a partial screening (half- hole). Indeed, 
the white line at the L3 edge is well reproduced in terms 
of shape and intensity, as well as the positions of the 
934, 936 and 943 eV structures. The only remaining 
source of disagreement with the experiment is related to 
the underestimation of the background between the two 
edges. 



A deeper analysis of the spectrum reveals that, as ex- 
pected, the transitions occur to the states formed by the 
hybridization between the e g of Cu and the p states of 
O. Although the s states of Cu and O lie in the same 
energy range as the empty 3d states, they are not seen 
at the Cu L2,3 edges. Indeed we find that transitions to 
states with s character yield a negligible contribution to 
the spectrum. 



C. S L 2 ,3 edges in molybdenite (2H-MoS 2 ) 

The 2H polytype of MoS 2 (2H-MoS 2 ) is a layered com- 
pound which crystallizes in a hexagonal structure. Each 
layer is composed of a M0S2 unit and there are two layers 
per cell (6 atoms). Layers are bound together by weak 
Van der Waals interactions. 

While isotropic in the case of cubic compounds as Cu 
and CU2O, the E1E1 x-ray absorption signal depends 
on the orientation of the polarization in the case of 2H- 
M0S2. We calculated the two independent XANES spec- 
tra corresponding to incident light polarized along the 
principal axis directions, in the plane and parallel to the 
c axis, at the 1^2,3 edges of S. Once again, two distinct 
sets of calculations have been performed, with and with- 
out core-hole. 

The comparison with experimental data^ for the two 
orientations of the incident beam polarization is shown 
in Fig. [51 The experimental spectrum at £ \ \ c was de- 
duced from the two orientations published in Ref. l58l . 
Calculations including the core-hole are in better agree- 
ment with the experiment than the ones without hole, 
at both orientations. The agreement with the experi- 
ment is good, with the position and spectral width of 
peaks being reproduced at least at the main edge. We 
note nevertheless that the energy position of the pre-edge 
peak is slightly overestimated in the theoretical spectra 
by 1 to 1.5 eV. Second, the intensity of the pre-edge peak 
is underestimated by the calculation. We believe that 
both disagreements are a consequence of the inter-edge 
mixing^ due to the fact the S 1^3 edges are very close in 
energy (1.1 eV). Third, the intensity of the spectral fea- 
tures above 180 eV is significantly overestimated by the 
calculations, at both orientations, albeit the energy de- 
pendent convolution. We have currently no explanation 
for this disagreement. 

In order to reproduce the quantitative XNLD at the S 
L 2 , 3 edges in 2H-MoS 2 , effects of both 2p -> s (Al = -1) 
and 2p — » d (Al = +1) transitions need to be included 
in the evaluation of the cross section (Fig. [5]). Fig. [5] 
shows the decomposition of the 1^,3 absorption into the 
two angular momenta channels, for the two polarization 
directions. Contrary to the common assumption that at 
L2,3 edges transitions to s states are negligible, Fig. [5] 
proves they have a significant weight in a wide range of 
energies, albeit lesser than the one of the 2p — > d transi- 
tions. This behaviour is a consequence of the significant 
El matrix element between the core 2p and valence s 



states. We expect it to be a general feature of L 2 ,3 edges 
of elements belonging to the third period (Al to Ar), as 
suggested from weighted s and d DOS calculations com- 
pared with Si L 2 , 3 edges in ELNES/XANES spectra of 
silicon^ and silicates^, as well as from calculations of 
ELNES at the L 2j 3 edges of Ga in GaN^ In particular, 
transitions to s states have been reported previously for 
XANES at the S L3 edge in In 2 S3^ based on a analysis 
of DOS. 

While the difference between the two sets of curves 
featured in Fig. |H] is due to transitions to s states, one 
can see that the aforementioned difference depends on the 
polarization direction albeit the spherical character of the 
s states. This apparent inconsistency can be understood 
by noting that, while the A^ = — 1 transitions are indeed 
isotropic (Fig. [9]), the total spectrum also contains cross 
(with respect to the selection rule), non- negligible d - 
s terms which do have an angular dependence. From 
Figs. [5] and [5] it is clear that the latter are an essential 
ingredient for describing the XNLD: in the absence of the 
cross term, the angular dependence of the 173 and 176 
eV structures would nearly disappear. 

We draw the attention on a fundamental difference be- 
tween K and higher rank edges. At K edges, the inten- 
sity of XANES structures is directly proportional to the 
corresponding I projection of the local DOS seen in the 
spectroscopy - the p-DOS (d-DOS) in the case of E1E1 
(E2E2) transitions. This is no longer the case at L2,3 
edges, in the presence of d - s cross terms, when the spec- 
tra can no longer be interpreted in terms of pure 1-DOS. 
Note that in the isotropic limit, the d-s interference term 
is strictly zero^ In this sense, the methodology of inter- 
preting polarized XANES spectra in non-cubic samples 
by confronting them to weighted DOS is not justified. 

In this respect, a challenge for the calculation is to dis- 
entangle a particular aspect of the spectroscopy, i.e. the 
interference between the Al = 1 and A^ = — 1 channels, 
from the electronic properties of the material itself. One 
can recover information about the 1-DOS from the spec- 
troscopy calculation and group theory. To begin with, in 
the following argument we consider for simplicity that (i) 
there is no spin-orbit coupling on the 2p state and (ii) the 
local point symmetry at the absorbing S site is described 
by an abelian group, with real characters. 

Defining (tpf\ = (s\ + (d\, and \ipi) = \pi) where (s\ 
and (d\ are the s and d components of the empty va- 
lence states (p being forbidden by the El selection rule) 
and \pi) with i = x,y, z is the core state, the El matrix 
clement in Eq. [T]is: 

\(il>f\fim\ 2 = \{ S \V\p^ + \(d\f>\p l )\ 2 +2(d\V\p l )(p. l \T>\s) 

(10) 

We focus on the last term in Eq. [TU] that represents 
the d-s interference. It is worthwhile to recall that a 
term of the form: 



{d\r a \pi)(pi\r a \s) 
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FIG. 8. (color online), x-ray linear dichroism at S L<2,3 edges 
in 2H-M0S2: calculations with (thick black solid) and without 
(red solid) core-hole versus experimental data from Ref. l5a 
(black symbols). The orientation of the polarization £ with 
respect to the c-axis is indicated in each panel. Calculations 
with core-hole were shifted with 1.5 eV to higher energies, 
with respect to those without hole. Calculations with and 
without hole were scaled independently to match the experi- 
mental data. 



FIG. 9. (color online) Decomposition of calculated XANES 
at the S L2,3 edges in 2H-M0S2 into angular channels. The 
total spectrum with core-hole (thick solid) is decomposed into 
transitions with Al = +1 (thin solid) and Al = — 1 (dashes). 
The orientation of the polarization £ with respect to the c-axis 
is indicated in each panel. The Al = — 1 channel is isotropic. 
The Al — +1 and Al = — 1 contributions do not sum up to 
the total spectrum, the remainder being due to cross terms. 



with a = x,y,z is zero under the above assumptions 
unless (i) a = i and (ii) d and s belong to the same irre- 
ducible representation. This can be understood by noting 
that, pi and r^ being in the same (one-dimensional) irre- 
ducible representation 7, 7 X7 is the totally symmetric ir- 
reducible representation A\. The term in pip is non-zero 
if and only if it contains the irreducible representation A\ , 
which only occurs (under the above assumptions) if (i) 
and (ii) hold. 

The point group of the S site in 2H-M0S2 is C^. In 
its abclian subgroup C s , s, on one hand, and d z i, d x 2_ y 2 
and d xy orbitals on the other hand, belong to the same 
irreducible representation and, as argued in the previous 



paragraph, can mix. One can easily see that in the £ \ \ c 
direction the cross term is d z i - s, while for the £ _L c 
orientation the cross term is due to d x i_ y i - s and d xy 
- s contributions. This conclusions remains valid if one 
considers the spin-orbit coupling on the 2p level. Under 
these circumstances, while the d yz - s and d xz - s contri- 
butions are not explicitly forbidden by symmetry, their 
contribution cancels after summation on the initial states 
belonging to the 2pi/ 2 or 2p\/2 manifold. 

The difference between the d - s interference in the 
two directions (the linear dichroism) is related to the 
anisotropy of the d orbitals. More specifically, the lin- 
ear dichroism of the d - s channel probes the difference 
between the d z 2 - s term, on one hand, and the joint con- 
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FIG. 10. (color online) Assignment of the main XNLD struc- 
tures to the d - s cross terms. Top panel: the cross terms for 
the two orientations, along c (thick solid) and in the plane 
(thin solid). Middle panel: the linear dichroism of the d - 
s cross terms, taken as the difference of the curves above. 
Bottom panel: calculated XNLD spectrum, obtained by the 
subtraction of the spectra corresponding to the polarization 
along c and perpendicular to the c axis respectively. The 
three panels share the same scale. The comparison between 
the middle and bottom panels proves that the XNLD mea- 
sures the anisotropy of the d states, via the mixing. The 
dichroic signal at 173 eV is mainly due to the d z i - s channel, 
whereas the structure at 178 eV comes from the d xy - s and 
d„i_,.i - s terms. 



tribution of d n .i_ 



s and gL 



on the other hand. 



This conclusion is equally valid had we considered the 
true C 3 i, symmetry. 

We have already seen in Fig. |H] that the main dichroic 
effects at 173 and 176 eV are due to the angular depen- 
dence of the cross terms. In Fig. [I0]wc continue the anal- 
ysis by comparing the calculated XNLD (bottom panel) , 
taken as the difference between the i \\ c and £ _L c spec- 
tra, to the cross d - s term in the two orientations (top 
panel) and its linear dichroism (middle panel) . The three 
insets sharing the same scale, it becomes clear that the 
XNLD can be safely interpreted in terms of d - s interfer- 
ence. Surprisingly, the pure AZ = 1 channel has no clear 
signature in the XNLD, except around 183 eV. While we 
have shown that the linear dichroism of the cross terms 
measures the d states anisotropy in 2H-M0S2, according 
to Fig. [TU] (middle and bottom panels) this is also the 
case for the XNLD spectra. One can therefore assign the 
dichroic feature at 173 eV to d z -i - s interference, while the 
feature at 176 eV is attributed to d x i_ y i - s and d xy - s 
mixing. In spite of the peculiar interference effect, XNLD 
still measures in-plane / out-of-planc d anisotropy. 

On the basis of the angular dependence, the authors of 
Ref. [58| assign the pre-edge peak and the shoulder at 170 
eV to transitions to s states. The decomposition in Fig. 
IH1 contradicts this conclusion: while there exists a con- 
tribution of the s states, it is definitely not dominating. 



A careful analysis of the DOS (not shown) shows that 
at the energies corresponding to the pre-edge peak the 
s states of S hybridize with the d orbitals of Mo. Alto- 
gether this proves the importance of calculations in order 
to accurately interpret the spectroscopic data. 

We recall that the interpretation of L 2j3 x-ray absorp- 
tion spectra in terms of single particle s and d orbitals is 
not exact in the presence of inter-edge mixing. However, 
we do not expect it to dramatically affect the main edge, 
where the probed d states are rather delocalized. This, 
together with the fact that the agreement with the ex- 
perimental data is satisfactory, re-confirms the choice of 
our methodology. 

The d - s cross term has also been depicted at the 
^2,3 edges of Ag in magnetic multilayers^ 3 - albeit not 
in the XANES spectra but in the x-ray magnetic circular 
dichroism (XMCD). This suggests that a correct descrip- 
tion of the cross term is equally important for magnetic 
spectroscopy. The spectroscopic interpretation would be 
even richer in the magnetic case, due to the exchange 
splitting within the 2pj / 2 and 2p 3 / 2 manifold. 



V. CONCLUSIONS 

In this work we have developed a technique based on 
DFT and on the PAW formalism in order to obtain the 
x-ray absorption cross section at a general edge, in both 
El and E2 approximations. We applied the method to 
the calculation of Cu L 2j3 edges in fee Cu and Cu 2 0, 
as well as to S L 2)3 edges in 2H-MoS 2 . In metallic Cu 
we find a good agreement with experimental data with- 
out having to include a core-hole. On the contrary core- 
hole effects are relevant in insulating Cu 2 and 2H-MoS 2 , 
where their inclusion is essential in order to obtain a good 
description of the experimental data. 

In the case of Cu L 23 edges in Cu 2 we equally study 
the effect of the Hubbard U term on the x-ray absorption 
spectra. In contrast to results on other oxides, we find 
that inclusion of U worsens the agreement with experi- 
mental data. 

Finally in the case of S L 2j3 edges in 2H-MoS 2 , we show 
that transitions to s states yield a non-negligible contri- 
bution to the S L 2]3 spectra in 2H-MoS 2 , in a wide range 
of energies. We believe that this is a general property 
of L 2:3 edges of third row elements from Al to Ar. We 
equally point out that when cross d - s terms are signifi- 
cant, the L 23 spectra are no longer to be interpreted in 
terms of projections of the DOS. The mixing is essential 
as to interpret the XNLD at S L 23 edges in 2H-MoS 2 . In 
such case it is utmost important to disentangle the two 
channels, which can only be achieved via first principles 
calculations. 

We have proven that single-particle approaches can 
be very useful for the interpretation of L 23 edges when 
the inter-edge mixing is negligible. In this case ab initio 
XANES calculation can assign spectral structures to var- 
ious angular moment channels, and explore their possible 
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interference. 



Appendix A: The El and E2 transition matrix 
elements 



For completeness we provide the full analytical expres- 
sions of the transition matrix element (<fi~R P \0\ipi) , at a 
general edge, as it was implemented in the XSpectra 
module of the Quantum Espresso distribution. 

By considering the spin-orbit coupling, the initial state 
(the core- level) \i[>i) is written as: 



i'n.hjmj (*") = ^ ^ (km. l , S1Tl s | jlUj 

m, = -!i m B = — s 



with the El transition matrix element becoming: 



(Al) 



(</>Rp|£> IVW.jm,) 



-in 



h s 

E E* 

li--^ ni s — — s 



m,j ,mi-\-m s \H''*'J 



{liirij — a, sa\jrrij} 



x E (-!)"£ E Y i"^)x s ms x s : 

fi— — 1 <r lm 

x [ dfY^(f) Yf(f) Y t m *(f) 
x dr r 3 R ni iimA r ) R P ia(r) 



By using the orthogonality of the spin functions: 

EE xtn.xy = EE <w 

m s a rn s a 



(A5) 



(A6) 



and the expression for the Gaunt coefficients (the integral 
is over the solid angle): 



where n^, li, nii, s — \/2 and m s are the usual nota- 
tions for the quantum numbers associated to the initial 
(core) state in the uncoupled basis and {hrrii, sm s \jm,j) 
the Clebsch-Gordan coefficients, whose tabulated values 
can be found in Ref. [6J. The functions R ni iim s { r ) an d 
Y ; m (f) are the radial and angular (complex spherical har- 
monics) parts, respectively, of the atomic wavefunction. 
Xm is the m s component of the spinor of spin s. The 
total momentum quantum numbers j = l{ — s, ...li + s and 
rrij = —j, ...,j define the core- level in the I ■ s coupled ba- 
sis. The Clebsch-Gordan coefficients are zero unless the 
following holds: 



rrij = m s + rrii 



(A2) 



The all-electron partial waves 0R P (r) are chosen as so- 
lutions of the Schrodinger equation of the isolated atom. 
As such they are written as 



Mr)=E E RMr)Yr{f)x S o 



(A3) 



where p labels the radial wavefunction at a given energy 
and n, /, m, s, a are angular and spin quantum numbers 
associated to the partial wave. 

The electric dipole transition operator can be ex- 
panded in the spherical harmonics basis as 



V = e-r = ^-r V (-1)" Y( l (f)Y-^(e) (A4) 
3 a— i 



0WM= / drYr*(r)Y™>(r)Yf(r) 



(A7) 



the El transition matrix element finally gives: 

(0R.pl V H'ndijnij) = y E ( ltm 3 ~ ° ' S<J ^ m o) 

a 
1 

x £ (-iy ¥{■>& e s'rv^ 

/j— — 1 lm 

x f drr 3 R nzha {r) R p i a {r) (A8) 

Jn R 

where f2# is the augmentation region. 

The electric quadrupole transition operator can be 
written as: 



Q = -(e-r)(k-r) = 

= l fc (T r V E W)lT(S)£ Y i X *(^ Y i X Ck) 



2 V 3 



/x=-i 



A=-l 



(A9) 



so that the E2 transition matrix element is 

(4>Rp\ Q \lpnilijmj) = 



2 V 3 



p; A 



x E ( liTn i ~ ^MJmj) Y^{f) G{^ x 



li rrtj~a,lfi 



ma 
R 



x / dr r 4 Rn iUrT (r) Rp la (r) (A10) 

Jo 
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The generalized Gaunt coefficients arc defined as: 



G 



2'< r n2- l lAmA 

'imi,l3m 3 



dO l™*(fi) Y™(fi) Y^*(Q) Y™(n) 



(All) 



and the following relation holds: 






Note that to calculate the cross section in Eq. [2] the 
contribution from each initial state ifii from the 2p\/2 
(2p 3 / 2 ) must be evaluated independently. Initial states 
with distinct jrrij quantum numbers do not cross in the 
single particle picture. 
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